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Abstract 

The symmetry and resonance properties of the Fermi Pasta Ulam chain with periodic 
boundary conditions are exploited to construct a near-identity transformation bring- 
ing this Hamiltonian system into a particularly simple form. This 'Birkhoff-Gustavson 
normal form' retains the symmetries of the original system and we show that in most 
cases this allows us to view the periodic FPU Hamiltonian as a perturbation of a nonde- 
generate Liouville integrable Hamiltonian. According to the KAM theorem this proves 
the existence of many invariant tori on which motion is quasiperiodic. Experiments 
confirm this qualitative behaviour. We note that one can not expect it in lower-order 
resonant Hamiltonian systems. So the FPU chain is an exception and its special fea- 
tures are caused by a combination of special resonances and symmetries. 

Keywords: periodic FPU chain, symmetry, resonance, Birkhoff-Gustavson normal 
form, near-integrability, KAM theorem 



1 Introduction 

The n particles FPU chain with periodic boundary conditions is a model for point 
masses moving on a circle with nonlinear forces acting between the nearest neighbours. 
It is in fact the n degrees of freedom Hamiltonian system on M 2n induced by the 
real-analytic Hamiltonian 



in which V : M — * M is a real-analytic potential energy function of the form 

V(x) = ^ + ^ + ^ + ... . (1.2) 
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The a, (3, . . . are real parameters measuring the nonlinearity in the forces between the 
particles in the chain. 



Numerically, the FPU system was first studied by E. Fermi, J. Pasta and S. Ulam, 
see These authors used the chain as a model for a string of which the elements 
interact in a nonlinear way. They expected that in the presence of small nonlinearities, 
the chain would show ergodic behaviour, meaning that almost all orbits densely fill up 
an energy-level set of the Hamiltonian. Ergodicity would eventually lead to an equal 
distribution of energy between the various Fourier modes of the system, a concept 
called thermalisation. FPU's nowadays famous numerical experiment was intended to 
investigate at what timescale thermalisation would take place. The result was aston- 
ishing: it turned out that there was no sign of thermalisation at all. Putting initially 
all the energy in one Fourier mode, they observed that this energy was shared by only 
a few other modes, the remaining modes were hardly excited. Additionally, within 
a not too long time the system returned close to its initial state. On increasing the 
strength of the nonlinearity, this recurrence occurred even earlier. Later computations, 
e.g. described in [0], confirmed that the same phenomena can also be observed in very 
large periodic chains. Empirical evidence was found that for small total energy, normal 
mode energies are hardly shared. Ergodic behaviour can only be observed when the 
energy level passes a certain critical value. 

In 1965 an article of Zabuski and Kruskal appeared, cf. 17]. These authors consid- 
ered the Korteweg-de Vries equation as a continuum limit of the FPU chain and nu- 
merically found the first indications for the stable behaviour of solitary waves, thereby 
suggesting an explanation for the striking data of the FPU experiment. In 1967, Gard- 
ner, Greene, Kruskal and Miura (Q) discovered infinitely many conserved quantities 
for the KdV equation, which should account for the regular behaviour of its solutions. 
Reference ]l(| contains a good overview of these results. They are suggestive, but do 
not provide a full explanation of FPU's observations as the impact of the transition 
from a discrete to a continuous chain has never been analysed. 

There is another, possibly correct explanation for the quasiperiodic behaviour of 
the FPU system. It is based on the Kolmogorov-Arnol'd-Moser (KAM) theorem (cf. 
0) and different from the Zabuski-Kruskal argument, it should work especially well for 
chains with a low number of particles. As is well-known (cf. 0), the general solution 
of an n degrees of freedom Liouville integrable Hamiltonian system is constrained to 
move in an ra-dimensional torus and is not at all ergodic but periodic or quasiperiodic. 
The KAM theorem states that most of the invariant tori of a nondegenerate integrable 
system persist under small Hamiltonian perturbations. Thus many authors, starting 
with Izrailev and Chirikov in Q, have stated that the KAM theorem explains the ob- 
servations of the FPU experiment. This reasoning seems plausible, but, as was clearly 
pointed out by Ford in ||, it is still completely unclear why the FPU system should 
be a perturbation of such a nondegenerate integrable system. This gap in the theory 
was recently mentioned again in the book of Weissert (|i~6|l). 

What does 'nondegenerate' mean here? Let us consider the frequency map u, 
which assigns to each ra-dimensional invariant torus of a Liouville integrable system 
the ra-dimensional vector of frequencies of the (quasi)periodic motion on this torus. An 
integrable system is called 'nondegenerate' if u; is a local diffeomorphism. The KAM 
theorem holds for perturbations of these nondegenerate integrable systems. 
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But it is no exception for an integrable system to be degenerate. A common exam- 
ple is the harmonic oscillator of which the frequency map is constant: the harmonic 
oscillator is highly degenerate. And indeed, perturbations of it are known that are 
ergodic even on low-energy level sets of the Hamiltonian. Ford gives a nice example 
of such a perturbation in his review article ||. We conclude that, although the FPU 
Hamiltonian can be considered as a perturbation of an integrable system -namely the 
harmonic oscillator-, the KAM theorem does not apply here! 

The aim of this paper is to overcome this problem. The method we use to do so is 
called Birkhoff-Gustavson normalisation - it is sometimes also called resonant normal- 
isation. It provides a transformation of phase space that in many cases enables us to 
write the periodic FPU Hamiltonian as a perturbation of a nondegenerate integrable 
Hamiltonian. 

It must be stressed that it seems highly exceptional that one can do this for a res- 
onant Hamiltonian system such as the periodic FPU chain. The current paper intends 
to make clear that the special symmetry, eigenvalue and resonance characteristics of 
the periodic FPU system play a crucial role in the construction of the near-identity 
transformation. It turns out that these characteristics cause the nondegenerate near- 
integrability of the chain. The conclusion is that the KAM theorem applies because 
of these resonance and symmetry properties: the quasiperiodic behaviour that Fermi, 
Pasta and Ulam observed is in some sense an exceptional feature of the FPU system. 



1.1 Outline of the paper 



This paper is a continuation of |12|] in which normal forms of small chains are computed 
and the KAM theorem is verified. We generalize and explain the results of [12] in this 
paper. 



In sections [2||6| the necessary theory is formulated. We start with an investigation 
of the eigenvalues (section and the discrete symmetries (section |||) of the periodic 
FPU chain. The concept of a Birkhoff-Gustavson normal form as an approximation of 
a Hamiltonian system is explained in section [|. It will be shown that normal forms for 
the periodic FPU chain exist that inherit its symmetry properties. 

In the appendix, which is based on notes of Beukers, number theory is used to 
compute all lower order resonances in the eigenvalues. We exploit this in sections [7| 



and g to prove theorem 8.2, which forms the core of this paper: it gives the restrictions 
that the Birkhoff-Gustavson normal form of any Hamiltonian with the same eigenvalues 
and symmetries as the periodic FPU chain, is subject to. 

These restrictions on the normal form allow us to point out many near-integrals 
of the chain in section |9[ We finish with an analysis of the /3-chain, which is proved 
to be near-integrable in section |n]. The KAM nondegeneracy condition can easily be 
checked when the /3-chain contains an odd number of particles. Some open questions 
are formulated for the even /3-chain. 
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2 Phonons 



To establish the sign conventions that we shall stick to during our analysis, some basic 
definitions follow here. For further reading on Hamiltonian systems and a thorough 
explanation of these concepts, the reader is referred to fl]]. 

We shall be considering Hamiltonian systems of differential equations on R 2n , the 
elements of which are denoted by (q,p) = • • ■ , q n ,Pi, ■ ■ ■ ,Pn)- On R 2n the sym- 
plectic form a := Yl]j=i dQj A dpj is defined. Endowed with this symplectic form R 2n 
is a symplectic space. Any Hamiltonian function H : R 2n — > R induces a Hamiltonian 
vector field Xh on R 2n which is defined by <j(Xh, •) = dH. Furthermore, for any two 
Hamiltonians F and G the Poisson brackets are defined as {F, G} := o~(Xp, Xq) = 
dF-X G = -dG ■ X F . 

Keeping these definitions in mind, we now start our analysis of the periodic FPU 
chain: 



In order to facilitate the equations of motion induced by the periodic FPU Hamiltonian 
( |l . l| ) , we apply a well-known Fourier transformation (q,p) i— > (q,p). For 1 < j < 8 
define 

« = V"L cos(— )% , ^ = A/ - 2^ «*(—)?* , 
fc=i '' fe=i 

. ,2jfcvr _ /2f . .2jfc7r 

^-i = \/-L sm(-^-)?* . Pn-i = v - 2^ sm (-^-)p* • 

Furthermore, define 



n ' — ' n V n * — ' n 

fc=l fe=i 



2.2) 



and if n is even, 



fc=i v k=l 

The new coordinates (q,p) are known as 'phonons'. The transformation to phonons 
is symplectic, that is a = Y^j=i^j A dp~j- For a proof, cf. |ll]] or [12]. In phonon- 
coordinates, the Hamiltonian reads 
n 1 

^"Eo^ 4 " + F 3(gi, • • • , 9n-i) + Hi(qx, ... , g n _i) + . . . , 

j=i (2.4) 

in which i/^ (fc = 2,3, ... ) denotes the fc-th order part of H; for j = 1, . . . , n, the 
numbers Uj are the eigenvalues of the linear periodic FPU problem: 

ua := 2sin(— ) . (2.5) 
n 

Exact expressions for H3 and in terms of the q~j can be found in the literature, cf. 
1 11]. We do not repeat them. 

The linearised equations are the equations induced by Hi- They read: 

Qj ' = Pj ,Pj ' = -Wj<ij ■ (2-6) 
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The Qj,pj (1 < j < n — 1) are harmonics with frequency loj; p n is constant, whereas 
q n increases with constant speed -note that oj n = 0. In fact, the linearised equations 
are Liouville integrable, with integrals Ej := \{p 2 + oj 2 q 2 ). The nonlinear equations (a 
or (3 unequal to zero) are much harder to analyse. The Ej are for instance no longer 
constants of motion. 



3 Reduction of a continuous symmetry group 



From (2.4) we see that H is independent of q n , even if a, f3, . . . ^ 0. This implies that 



p n is an integral of H. The set ^^{O}) defines a 2n — 1 dimensional hyperplane in 
M 2n , invariant under the flow of both Xu and Xp n . The flow of Xp n = induces a 
symplectic M-action on this hyperplane. The time-t flow e txpn is actually given by 

e "" n - §(*4 +ft 4) +t i ' (3 - 1) 

or written out in the original coordinates: 



The orbits of this flow are the lines (q,p) H-K^p-. It is clear that the 2n — 2 dimensional 
hyperplane g^ 1 ({0})np~ 1 ({0}) = M 2n ~ 2 is transversal to these orbits. Therefore, R 2n ~ 2 
is a model for the space P^" 1 ({0})/k of Ap n -orbits lying in p~ 1 ({0}). R 2n ~ 2 inherits the 
symplectic structure a := X]j=i ^Qj ^ dpj from M? n . And since the FPU Hamiltonian 
H is constant on the orbits of the flow of Xp n , H reduces to a Hamiltonian on M 2n_2 
given by 

n— 1 ^ 

# = 9 ^ + + ^fe' • • • ' ^n-l) + • • • , 9n-l) + • • • • 

j=i (3.3) 



The reduced Hamiltonian (|3.3| ) represents the periodic FPU system from which the 
centre of mass motion has been eliminated. 

Since uj 2 >0(l<j<n — 1), we conclude with the Morse-Lemma (cf. [jlj) that 
the level sets of H are 2n — 3 dimensional spheres around the origin of M 2n_2 . And 
since H is a constant of motion for the flow of Xh, we see that the origin is a stable 
stationary point for the reduced system induced by the reduced Hamiltonian (3.3). 



4 Discrete symmetries 



Apart from the continuous family of symmetries of the previous section, the FPU 
Hamiltonian has some discrete symmetries. These have important dynamical conse- 
quences. 
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The first discrete symmetry is a rotation symmetry. Let T : R 2n — ► R 2n denote the 
circle permutation, the unique linear mapping defined by 

m d d d d , A . 

T : i — > , . (4.1) 

dqj dqj-x dpj dpj-i 

T is symplectic: T*a = a. Furthermore, note that T leaves H invariant: T*H := 
H o T = H. This implies that the Hamiltonian vector field Xh induced by H is 
equivariant under T : DT ■ Xh = Xh ° T. In other words: if 7 : R -> R 2n is an 
integral curve of Xh, then T o 7 : R — ► R 2n is an integral curve of Xh- This is why 
we call T a symmetry of H . The same thing holds for the powers of T. The group 
(T) := {Id, T, T 2 , . . . , T n ~ 1 } = Z/ n z is a discrete symmetry group of H. 

We can point out another discrete symmetry, namely the reflection S : 
which is the unique linear mapping sending 



[>2n . in>2n 



d d d d 

S : — h+ — , — ^ — . (4.2) 

dqj dq n -j dpj dp n -j 

S is again a symplectic symmetry: S*a = a and S*H = H. The group (S) := {Id, S} 
= Z/2Z, whereas the full discrete symmetry group {T,S) := {Id, T,T 2 , ... ,T n , 
S,ST 2 , . . . , ST™" 1 } = D n is called the 'n-th dihedral group'; its group structure is 
determined by the relation ST = T n ~ 1 S. The vector field Xh is equivariant under the 
elements of (T, S), that is (T, S) maps integral curves of Xh to integral curves of Xh- 

The reader should note that T and S leave 9^ 1 ({0}) Hp~ 1 ({0}) invariant. There- 
fore, T and S reduce to linear symplectic mappings on R 2n_2 that leave the reduced 
Hamiltonian invariant^]. 



5 Normalisation 

We shall study the reduced FPU system ( |3.3[ ) using Birkhoff-Gustavson normalisa- 
tion. In fact, we shall construct a near-identity transformation of phase-space allowing 
us to write the FPU Hamiltonian in 'normal form', meaning that it can be seen as 
a perturbation of a rather simple system. The study of the truncated normal form 
-that is this simpler system- leads to important conclusions for the original FPU sys- 
tem. For instance, the solutions of the truncated normal form are approximations of 
low-energetic solutions of the original system valid on a long time-scale. Integrals of 
the truncated normal form are near-integrals of the original system: on orbits of low 



1 The FPU Hamiltonian also has a reversing symmetry, namely the mapping R : K 2n — > K 2 ™ given by 

dqj dqj ' dpj dpj 

R leaves the FPU Hamiltonian invariant, i.e. R*H = H. R is anti-symplectic in the sense that R*a = —a. 
This implies that the vector field Xh is anti-equivariant under R: DR ■ Xh = — Xh R- In other words: if 
7 : K — > R 2n is an integral curve of Xh, then R o 7 o (—Id) : K — > R 2 ™ is an integral curve of Xh- Since R 
leaves g^" 1 ({0}) flp~ 1 ({0}) invariant, R reduces to an anti-symplectic mapping on R 2n ~ 2 leaving the reduced 
Hamiltonian invariant. More information on reversing symmetries can be found in B. 
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energy, they are almost conserved for a long time. See []T4| for an explanation and ex- 
plicit statements. Furthermore, the truncated normal form can help us understanding 
bifurcation phenomena. And last but not least, if the truncated normal form of the 
FPU chain is integrable in a nondegenerate way, then the FPU chain is a perturbation 
of a nondegenerate integrable system. We may apply the KAM theorem then and 
conclude that almost all low-energetic solutions of Q3.3|) are quasiperiodic and move on 



tori. Conclusions of this type were drawn for the first time in [12] 



The setting of normalisation is the following: 

Let Pk be the set of all homogeneous fe-th degree polynomials in (qi, . . . , q n -i,pi, 
... ,p n -i). The set of all power series without linear part, P := © fc>2 -Pfc! is a 
Lie-algebra with the Poisson bracket. For each h G P the adjoint representation 
ad/j : P — ► P is the linear operator defined by &dh(H) = {h, H}. Note that whenever 
h G P k , then ad h : P t -> P k +i-2- 

The flow e tXh of a Hamiltonian vector field Xh induced by h G P — P2 is a symplec- 
tic near-identity transformation in M? n ~ 2 . For its action on an arbitrary Hamiltonian 
H G P we have ^(e tXh )*H = dH ■ X h = —adh(H). This is a linear differential equa- 
tion in P of which the solution is (e tXh )*H = e~ ta,dh H. In particular the near-identity 
'Lie-transformation' e~ Xfl = Id — Xh + ■ • • transforms H into 

H' := (e- Xh )*H = e adh H = H + {h,H} + ^{h,{h,H}} + ... . (5.1) 

Let us denote the A:-th order part of the Hamiltonian H -that is the projection of H on 
Pk- by Hk- If for instance h G P3, then we obtain the formula H' k = Y^n=Q ijjr (ad^)" 1 
(Hk- m ). We just gathered all terms of equal degree in formula (|5.1[). 



Assume now, as is the case for the reduced FPU Hamiltonian, that ad# 2 : Pk — > Pf. is 
semisimple (i.e. complex-diagonalisable) for every k > 2. Then P k = ker ad# 2 © 
imad#2- In particular H3 is uniquely decomposed as H3 = 73 + 53, with G 
ker ad# 2 , 53 G im ad# 2 . Now choose a /13 G P3 such that ad.H2i.h3) = 53. One could for 
example choose hs = g$ := (ad# 2 |- m ac j^ ) (sa). But clearly the choice /13 =53+^3 

suffices for any P3 G ker adn 2 H P3. For the new Hamiltonian H' we calculate from 
( |5.1| ) that H' 2 = H2, ^3 = /3 6 ker ad# 2 , H' 4 = H4 + {I13, H3 — ^53}, etc. But now we 
can again write H' A = f± + 54 with f± G ker ad# 2 > ff4 G im ad# 2 and it is clear that by 
a suitable choice of /14 G -P4 the Lie-transformation e _Xh 4 transforms our H' into ii"" 
for which H' 2 ' = H2, H' 3 ' = f$ G ker ad# 2 and H'l = f± G ker ad# 2 . Continuing in this 
way, we can for any finite r > 3 find a sequence of symplectic near-identity transforma- 
tions e~ Xfl 3 , . . . , e~ Xhr with the property that e~ Xfl k only changes the H\ with / > k, 
whereas the composition e~ Xhr o . . . o e~ Xh s transforms H into H with the property 
that Hk Poisson commutes with H2 for every 2 < k < r. H is called a normal form of 
of order r. Its study can give us useful information on low-energetic solutions of the 
original Hamiltonian H. More on normalisation by Lie-transformations can be found 
in 01. 
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6 Normal forms and discrete symmetry 



In section |] we investigated the discrete symmetries of the periodic FPU Hamiltonian. 
We saw that they reduce to symmetries of the reduced FPU system on M? n ~ 2 . In this 
section we show how one can construct normal forms of the reduced FPU Hamiltonian 
that have the same symmetry properties as the reduced FPU Hamiltonian itself. This 
will help us formulating some important restrictions on the normal form in section ||. 

The symmetry properties are captured in the definition of the symmetric subspace 
of P: 

P ST : = {/ G P\ S*f = f, T*f = /} . 

Note that the FPU Hamiltonian is in P . 

The next observation is that S* and T* are Lie-algebra automorphisms of P: 

S*{f, 9} = {5*/, S*g} , T*{f, g} = {T*f, T*g} . (6.1) 

simply because T and S are symplectic. Now take / G P ST and g G P ST . Then 
from Q it follows that S*{f,g} = {S*f,S*g} = {f,g} and T*{f,g} = {T*f,T*g} = 
{/, g}. This means that P ST is a Lie-subalgebra of P: if /, g G P ST , then {/, g} G P ST . 
Alternatively stated: if h G P ST , then ad^ : P ST — > P ST . In particular, e ac ^ : P ST — > 

pST 

Since ad# 2 leaves P ST invariant, we know that P ST = (ker ad# 2 n P ST ) © (im ad# 2 H 
P ). So if we decompose the third order part of the FPU Hamiltonian as H3 = fs + gs 
with /3 G ker ad# 2 , 53 G im ad# 2) then f^, g% G P^ T automatically. /13 = 53 = 
(ad^lixn ac l# ) _1 (fi l 3) is the unique element of im ad# 2 nP 3 ST for which ad# 2 (/i3) = 53. 

But since g 3 G P 3 5T , we find that H' = (e~ x ^)*H = e ad s 3 H G P ST . Of course the 
choice /13 = 53 +P3 also suffices for any ^3 G ker ad# 2 n P 3 ST . 

It should be clear that continuing this procedure, we can produce normal forms 
H G P ST of H up to any finite order^]. 



7 Simultaneous diagonalisation 

From ( |6.1|) we infer that 

(T* o adf/ 2 )(/) = T*{Pi2, /} = {T*H 2 ,T*f} = {H 2 ,T*f} = (ad H2 of)(/) . 

(7.1) 

So ad# 2 and T* commute on P^ . Therefore adn 2 leaves the eigenspaces of T* invariant 
and we can diagonalise ad# 2 and T* simultaneously. This allows us to calculate the 
subspace Pj. n ker adjy 2 n ker(T* — Id) C Pk in which is contained and helps us 
formulate some important restrictions on the normal form of the FPU Hamiltonian. 

In order to perform this simultaneous diagonalisation, we introduce the 'super- 
phonons' (z,Q. For 1 < j < §, define: 

2 Although the bookkeeping is a bit harder, one can extend the previous argument to prove that the 
normal forms can also be chosen invariant under R*. For a complete proof, cf. H. 



8 



1 • 1 n 

z i := o \Pi ~ % Vn-j) + -rriQj ~ Hn-j) = —t== V e ™ (pjfc + i^g*) 

2 2 V2nfrt 



Cj : = (Pj + iPn-j) - Tjilj + = : — ~f= X] e » (p k - iujjq k ) 



n-j 



2iujj 2 iujV2n— ^_ 2 ) 

~o(Pj ~ + - iQn-j) = 7= > e « (p fc - aw.,-?*.) 

2 2 ^ 2n 



1 1 1 " 2 fc 

Cn-i := w—,(Pj + iPn-j) + y(c[j + iq n -j) = -j= ^2 e ~^(Pk + iujq k ) 



and if n is even: 



1 1 

z ? := — + = - — 7^= y^(- i ) fc bfe + ?*) ( r - 3 ) 

2 2 K=l 

1 1 " 

Cf := ^|(pf " iwngn) = — Y^(-l) k (Vk - i^q k ) 

One checks that {zj, z k } = {Q, £&} = and {zj, ( k } = Sj k , the Kronecker delta. So our 
superphonons define canonical coordinates, i.e. a = Y^j=i z j ^ Cj- 



From ( [4.1[ ) we infer that T*qj = qj+i and T*pj = Pj+i, where qj,Pj '■ M. — > M 
are the coordinate functions. So from (7.2) we see that 

^ 2-Kij 2tv ij 2irij 2irij 

T : Zjh^e~Zj, Cj ^ e n Cj, z n -j^e~z n -j, Cn-j ^ e » Cn-j , 

Z2 i — > —zn and Cn | — > — • (7.4) 

22 2 2 vy 

We conclude that T* acts diagonally on (z, £)-coordinates. And it acts diagonally on 
monomials in (z, Q: if 0, 6 £ {0, 1,2,... } n_1 are multi-indices, then 

T* : A* - e 2 ^ 1 A* , (7.5) 

fx being defined as: 



M(e,0):= £ j(e i + e n ^-0 i -0„_ i ) + £(e»-0») modn . 

i<i<f ■ (7.6) 

On the other hand one calculates: 



H 2 = ^ iuJ j( z j(j ~ Zn-j(n-j) +iw*2|C§ • (7.7) 

i<i<f 

So we also diagonalised ad# 2 with respect to monomials: 
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ad H2 : z e C >-> v(@,9)z e (° , (7.8) 

in which v is defined as 

u(@,6) := iuJ j(8j ~ O n -j ~ 9j + @n-j) + - 0» ) . (7.9) 

i<i<? 

Monomials z^C 6 * commuting with H2 -the ones for which u(Q,9) = 0- are called reso- 
nant monomials. They are particularly important because they cannot be normalised 
away. 



8 Restrictions for symmetric normal forms 

From section ^ we know that we can transform the periodic FPU Hamiltonian into a 
discrete symmetric normal form of any desired order. Suppose we did so up to order 
r. Then H k G P k n ker ad# 2 n ker(T* - Id) for any 2 < k < r. But since both T* 
and ad# 2 act diagonally in (z, ^-coordinates, we know that this H k must be a linear 
combination of monomials z (, e for which 

|G| + \B\ = k , n(Q, 9) = mod n and v(Q, 9) = . (8.1) 

Extra restrictions on H k , with which we shall deal later, arise from the fact that H k 
can be chosen in the even smaller set P ST []. But first we investigate which and 

9 satisfy ( |8.1| ). Because the ujj in ( |7.9|) are of the form 2isin(^), this is actually a 
number-theoretical question that we shall solve for |G| + \0\ = 2, 3, 4. 

The quadratic case - i.e. |0| + \6\ = 2 - is easy: since all the Uj are different, we 
find from v{®, 9) = that the Lie-subalgebra P2 Pi ker ad# 2 C P2 is spanned by the 
monomials 

ZjCj, Zn-jCn-j, ZjZn-j, (j(n-j (1 < j < ^) aIld Z f Cf ■ (8-2) 

T* acts diagonally on these basis-elements as follows: 

T* : ZjCj I ^ ZjCj, ^n-jCn-j ^ ^n-jCn-j, Z| Cf ^ Zf Cf , (8-3) 
z j z n—j 1 e n z j z n—ji CjCn~j 1 * e n CjCn~j ■ 

The Lie-subalgebra i-^Hker ad# 2 Piker (T*— Id) = spanjzjCj, z n-j(n-j, z f Cf } is abelian. 
From (|4.2| ) and (7.2) we calculate the action of S* on the coordinate-functions: 

S* : Zj (-> -lUJjCn-j, Cj ^ " Zn-j ^ iw jCj, Cn-j ^ — ^2 ^ ~ z %, C# ^ ~ 



3 and invariant nnder R* 
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So the action on the basis-elements reads: 



S* : ZjCj H-> -Z n -jCn-j, Z n -j( n -j l-> ~Zj(j, 0|C§ ^ Cf ) (8-5) 
ZjZn—j 1 > ^jCjCn—ji CjCn—j 1 2~ZjZ n — 3 ' 

We conclude that the Lie-subalgebra n ker ad# 2 is spanned by the quadratics 
ZjCj — Zn-jCn-j and Z|C 2 - Note that H2 itself is indeed a linear combination of these 
quadratics. 



The analysis is harder if we consider the cases |0| + \9\ = 3, 4. With the use of 
number theory, the proof of the following theorem is given in the appendix. 

Theorem 8.1 

i) The set of multi-indices (G, 9) € {0, 1,2,... } 2n ~ 2 for which |0| + \0\ = 3, /x(0, 9) = 
mod n and 9) = is empty. 



ii) The set of multi-indices (6, 9) £ {0, 1,2,... } 2n ~ 2 for which |6| + \9\ = 4, /x(0, 0) = 
mod n and z^(0, 9) = is contained in the set given by the relations 9j — 9 n -j — @j + 
Q n _ j = 0n - 0n = 0. 



Theorem 8.1 has some major implications. We shall investigate these now and they 



will be summarised in theorem 8.2 



From i) we see that P 3 5T n ker ad# 2 C P 3 n ker ad# 2 n ker(T* - Id) = {0}. 

First of all, this implies that we can always transform away H3 from the peri- 
odic FPU Hamiltonian: H% = 0. This is an unexpected result. Consider for ex- 
ample the chain with 6 particles, which satisfies a third order resonance relation: 
u>\ : u>3 '. c<->5 = 1:2:1. For systems with a third order resonance relation one 
can generally not expect H% to be trivial. But, as was observed for the first time in 



12], it is trivial for the 6 particles chain. One could say that the 1:2: 1-resonance 
is not active at i^-level. We now know that for the periodic FPU chain no resonance 
will ever be active at H^-level. This simplification is caused by the symmetries of the 
FPU system. 

Secondly, we conclude from i) that the /13 of section |(| is uniquely determined by 
the requirement that it be in P^ T . This in turn uniquely determines H4. 



From ii) we infer that any element of P4 Pi ker ad# 2 n ker(T* — Id) must be a lin- 
ear combination of products of two of the basis-elements in ( |8.2|) . 

Note however that not all these products are really T* -invariant and that the full 
normal form is even invariant under S*. We work out these extra restrictions now. 

The question which products of the basis-elements (|8.2D are invariant under T* is 
easy to answer with help of the formulas (8.3). Clearly, all products of ZjCj> z n~jCn-j 
andz-Cf are. T* multiplies the terms (zjCj)(z k z n - k ), (^0)(CfcCn-ifc), (Zn-jCn-j)(zkZn-k), 
(z n -j(n-j)((k(n-k), (zn(n)(z k z n ^ k ) and (z^Cf )(CfeCn-fe) with a factor e ± ~^~ / 1, so 

4ni(j — k) 

these terms are not invariant under T*. T* multiplies fe^n— i)(CkCn— Jfc) by e ™ 
which is 1 if and only if 2(j — k) = mod n. But because 1 < j, k < S, the 
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condition is 2(j — k) = 0, i.e. j = k. Thus we end up with a term that we al- 
ready had: (zjZ n -j)(CjCn-j) = (zj(j)(Zn-jCn-j)- Finally, the terms (zjZ n -j)(z k z n - k ) 
and (CjCn-j)(CkCn-k) are multiplied by a factor e ™ which is 1 if and only if 
2(j + k) = mod n. But since 1 < j, k < ^, the only possibility is that 2(j + k) = n, 
that is n must be even and j + k = ^. This concludes our search for fourth order 
monomials invariant under T* and Poisson commuting with H2. 

We shall check now which combinations of these terms are also invariant under S* . 
The action of S* on P2Hker ad# 2 can be diagonalised in real coordinates. For this pur- 
pose, besides our familiar complex basis, we also define the following real basis-elements 
for P 2 n ker ad# 2 . For 1 < j < §, let 

■( A /• \ C-2 1 -2 1 2-2 , 2-2 \ 



Oj := i(2jCi + Zn-jCn~j) = PjQn-j ~ Pn-jQj , (8.6) 



i 1 

dj := —{UjCjCn-j ~ ZjZn-j) = —(pjpn-j +Ujqjq n -j) , 
UJj Uj 

and if n is even 

> / -2 2 -2\ 

2 2 S 2 2W» 2 2 2 

2 

Note that these basis-elements are subject to the relation 

a) = b) + cf + d] (8.7) 
and that H 2 can easily be expressed as 

H 2 = u 3 a i • ( 8 - 8 ) 

Our definitions diagonalise the action of S*: 

S* : cij 1 ^ aj, o» wa», 6j- 1— > — 6j, Cj 1— > Cj, 1— > — dj . (8-9) 

The products aja^, a?±aj and Oj&fc are invariant under S* and T*. The products ajb k 
and a^&fc are not invariant under S* , although they are under T*. It is left as an easy 



excercise for the reader to prove that the only configuration for other terms to appear 
is djdik-j — CjCr±_j. 

We summarize the results of this section in the following theorem: 



Theorem 8.2 Let H be the reduced periodic FPU Hamiltonian (3.Z ). There is a fourth 
order normal form H of H which is invariant under T* and S* For this normal 
form we have H3 = 0, whereas H4 is a linear combination of the fourth order terms 
ajOk, bjbk (1 < j, k < ^) and if n is even an.aj (1 < j < ^) and djdn.-j — CjCn_j 
(l<i< f)- 



'and R* 
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9 Near-integrals or integrals of the truncated 
normal form 



In the previous section we proved that the truncated fourth order normal form of the 



periodic FPU Hamiltonian is subject to many restrictions, as indicated in theorem S.2 



This enables us to point out some integrals for the truncated normal form. These are 
near-integrals of the periodic FPU chain: quantities that are nearly conserved by the 
flow of the orginal chain ( |3,3| ) for a long time, cf. [14|. 

In order to be able to compute these integrals, we first write down the commutation 
relations between the real basis-elements (8.6). They are given by 



{bj,Cj} = 2dj , {bj,dj} = —2cj , {cj,dj} = 2bj . (9.1) 

All the other Poisson brackets between basis-elements give 0. These relations lead to 
the following conclusions: 

9.1 The odd chain 

Corollary 9.1 If n is odd, then the truncated normal form H2 + of the periodic 
FPU chain is Liouville integrable with the quadratic integrals aj,bj (1 < j < ^i— )■ 

Proof: H2 is a linear combination of the quadratics aj and H4 is a linear combination 
of the fourth order terms a^a^ and bjb^. The aj and b^ Poisson commute with all these 
terms and with each other. □ 

It is well-known (cf. Q), that the integrals of a 2n — 2-dimensional Liouville inte- 
grable Hamiltonian system generally define n — 1-dimensional invariant tori. Let's see 
what these tori look like here and analyse the integral map F : IR 2 " --2 — > R n_1 that 
maps (q,p) 1— > (a, b): 

Proposition 9.2 

im F = {(a, 6) G M™" 1 ^ > 0, \bj\ < a,} . (9.2) 

For (a,b) G {im F) int = {(a,b) G M. n - l \aj > 0,\bj\ < aj}, F' 1 ^^)}) is a smooth 
n — 1-dimensional torus . 

Proof: Clearly, im aj = [0, 00). The level set of aj is, for aj > 0, the cartesian prod- 
uct of M 2n_6 and a 3-dimensional sphere in R with radius \j2aj . Let us consider bj 
restricted to the level set of aj. To compute its critical points, we use the method 
of Lagrange multipliers: (q,p) is a critical point iff there is a constant A such that 
Daj(q,p) = XDbj(q,p), that is (ujjqj,u)jq n -j, ±pj, ±-.p n -j) = H-Pn-j,Pj,q~n-j, -Qj)- 
From this we infer that A = ±1. For A = 1, we find p n -j = ~ u; jQj^ Pj = ^jQn-j- In 
these points we have bj = aj. A = — 1 gives p n -j = ^jQji Pj = —^JjQn-j, so bj = —aj. 
These are the extreme values of bj on the level set of aj, giving ( |9.2[) . We also learn 
from this that if aj > and \bj\ < aj, then Daj and Dbj are independent. So if 
(a, b) G (imF) mi , then all Daj and Dbk are independent on F _1 ({(a, b)}). According 
to a theorem of Arnol'd (cf. such a level set must be a torus. □ 
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In order to compute the flow on these tori, we make the explicit transformation to 
action-angle coordinates (q,p) i— > (a, b, <fi, t/j) as follows. Let arg : M 2 — {(0, 0)} — ► K/2ttZ 
be the argument function, arg : (r cos r sin $) i-> Then, for 1 < j < ^i, define 

bj ■= PjQn-j - Pn-jQj , (9.3) 

h '■= \^&{-Pn-j ~ UjqjiPj - Ujq n -j) + -axg(p n -j - Ujq~j,pj + Ujq n -j) > 

1 1 

ipj ■= 2 ar g(-P"-i - Ujq^Vj - ujjqn-j) - -arg(p n _j - Ujqj,pj + Ujq n -j) ■ 



Note that these are well-defined as long as (a, b) G (imF) m( . With the formula 

cdy—yd 



d arg(x,y) = x ^ + ^ x , one can verify that the (a,b,4>,ip) are canonical coordinates: 
a = Y^j=i d<jj h dpj = J2i<j<n—l dcfrj A da,j + dtpj A dbj 



The truncated normal form H2 + H4 is a function of the actions aj , 6j . Its induced 
equations of motion therefore read: 

aj = bj = , (9.4) 

• dH 4 , . ■ 9i^4 . . 

^■=0,, + — (a,6), ^ = — (a,6), 

which are very simple. In order to verify that that the truncated normal form H2 + 
is nondegenerate, we examine the frequency map uo which adds to each invariant torus 
the frequencies of the flow on it: 



/ , \ / dH± dH ^ dH^ dH^ \ 

to : (a,b) \-> \ui + — — (a, b), . . . ,Wn=i + (a, b), (a, 6), . . . , — (a, b) . 

y oai 2 dan^ obi obn^ J 

f\'2~TJ rj2Tf 

u is a local diffeomorphism iff both the constant derivative matrices da .g^ and Q b .g^ 
are invertible. We will explicitly check this for the odd /3-chain in section |ic|. 

The situation is more difficult in the case of 



9.2 The even chain 

Corollary 9.3 If n is even, then the truncated normal form Hi + H4 of the periodic 
FPU chain has the quadratic integrals aj (1 < j < ^) and bj — b^_j (1 < j < j). 

Proof: H2 is a linear combination of the quadratics aj (1 < j < ^), whereas H4 is a 
linear combination of the fourth order terms aja^ (1 < j,k < ^), bjb^ (1 < j < ^) and 
djdi±_j — CjCi±_j (1 < j < j). The aj clearly commute with all these terms. So do the 
terms bj — bii^j: {bj — bn^j, bk} = {bj — b^-j, a^} = {bj — bii^j, a«} = 0. But one also 
verifies from (|9.1| ) that \bj — biL_j,CjCi±_j — djdii^j} = c?±_j{bj,Cj} — Cj{b™__j,cii_j} 
-d^_j{bj,dj} + dj{bn_j,dn^j} = 2cn_jdj - 2 Cj dn_j +2d^_j C j - 2d i c«_ j = 0. □ 
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If n is even, then the n-l-degrees of freedom Hamiltonian H2+H4 has at least 3n ^ 4 (if 4 
divides n) or 3 " 4 ~ 2 (if 4 does not divide n) quadratic integrals. These are near-integrals 
for the original chain ( |3.3[ ). We have not yet found a complete system of integrals for 
the truncated normal form though. We will do so for the even /3-chain in section 10. 2j . 



10 The normal form of the /3-chain 

In this section we present the explicit normal form of the periodic FPU Hamiltonian 
in the case that H3 = 0, i.e. a = in (p..lj). This chain, that has no cubic terms, 
is usually referred to as the /3-chain. A calculation of the normal form of order 4 is 
relatively easy in this case, because one does not have to transform away #3 first. The 
calculation is still tedious though and that is why we do not present it. The reader can 
find an example of a similar computation in [|l2j |. The following theorem is a major 
generalisation of the result in |l2j , which in turn is a restatement -with a much more 
efficient proof- of a theorem in the PhD thesis of Sanders (|i~3|1). 

Theorem 10.1 If a = 0, then in the periodic FPU chain one has 

— /3 v UJkUJl , W L Q 2 12 \ , 1 2 , 1 \ ^ 

Hi= n\ ^ — a kai+ 2^ ^( 3 «fc- 6 fc) + 4«f + 2 a f Z> Wfcflfc 

[o<fc<Z<§ 0<fc<§ 0<fc<§ 

8 <i^ 16 4 4 J (10.1) 

In formula (10.1) it is understood that terms with the subscript | and § only appear 
if 2 respectively 4 divides n. 



10.1 The odd /3-chain 

In formula (10.1) we see again what was already predicted in theorem |8.2| , namely that 
H4 is a linear combination of the terms a^a k and bjb k (1 < j,k < ^j— ). According to 
corollary |9.1| this normal form is integrable, the aj and bj being the (quadratic) inte- 
grals. To check the nondegeneracy condition, we compute the second order derivative 
matrices of H4 with respect to the action variables aj and bj: 



d 2 H 4 



p 



( M 

4lx>2 0J\ 



dajda k 16 



n 



AUJ1UJ2 



\ 2 2 



2 
2 



3w|_i 1 
2 / 



(10.2) 



d 2 H 4 
dbjdb k 



JL 

16 n 



re-1 / 
2 ' 



(10.3) 
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g b ^ k is clearly nondegenerate. But so is ® a ^ k ■ This can be proved by applying 
elementary row and column operations to ( [10. 2 ), thus reducing it to upperdiagonal 



form. This yields an expression for the determinant that is unequal to 0. We conclude 
that the reduced periodic /3-chain with an odd number of particles can, after a near- 
identity transformation, be written as a perturbation of a nondegenerate integrable 
Hamiltonian system. Therefore, the KAM theorem (cf. [§]) applies: 

Theorem 10.2 If n is odd, a = and (3^0, then almost all low-energy solutions 



of the reduced periodic FPU chain (3.8) are periodic or quasiperiodic and move on 
invariant tori. In fact, the relative measure of all these tori lying inside the small ball 
{0 < H < e], goes to 1 as e goes to 0. 

It should also be possible to write down an expression for the normal form if a ^ 0. 
The nondegeneracy condition can be checked again then. But the computation of 
this normal form is very nasty - transforming away H3 we obtain the transformed 
H' 4 = H4 + 7j{(ad# 2 |j m a( ^ H )^ 1 (H^), H^} which thereafter has to be normalised to 

produce H^. The result is most likely that for almost all a and (3 the nondegeneracy 
condition holds and the KAM theorem applies. Without computation this is clear for 
|a| <C because in this situation the coefficients of the normal form (10.1) change 
only slightly and the invertible matrices form an open set in the set of all matrices. 

10.2 The even /3-chain 

It is a surprise that in the normal form of the even /3-chain no terms bjbk (j 7^ k) arise, 
see formula (10.1). This leads to the following remarkable conclusion: 

Corollary 10.3 // n is even and a = 0, (5 7^ 0, then the truncated normal form 



H 2 + H4 of the reduced periodic FPU chain (3.5) is Liouville integrable. The integrals 
are the quadratics chj (1 < j < |), bj — 6»-j (1 < j < f ) and d™ (if n is a multiple of 

4) and the fourth order terms u\b\ + w\_ k b\_ k + 4w| fc (cfcC«_fc — d^d^-^) (1 < j < f ). 

22 22 

Proof: This follows from simply computing all the Poisson brackets, using ( |9.1[ ) and 
the fact that the Poisson brackets form a derivation. □ 

Only the a,-,6,- — bn_ i and dn_ induce a 2-7r-periodic flow and can therefore be seen 
as actions after some symplectic action-angle transformation. It is an open problem to 
construct the remaining action variables. Thereafter one could differentiate H4 with 
respect to them and verify the KAM nondegeneracy condition. 

One exceptional case is easier: the /3-problem with 4 particles. Its truncated normal 
form reads: 

H 2 + H 4 = V2a x + 2a 2 + j (^4 + X -a\ + ^-a x a 2 + j , (10.4) 
which has the commuting integrals ai,a 2 and d\. The frequency map is 



uj : (01, a 2 , di) h-> (V2 + — a\ H —a 2 , 2 + -a 2 H — ai, — di, 

lb 8 8 8 lb 



(10.5) 
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uj is nondegenerate, since 



d(ai,a 2 ,di) 




(10.6) 



is invertible. So a similar theorem as 10,2 holds for the /3-chain with 4 particles 



It is unclear what happens for the even chain if a ^ 0. The truncated normal form 
might not be integrable. On the other hand we already found about ^ integrals. 
And in [|12| it was already shown that the normal forms of the a-/3-chain with up to 6 
particles are Liouville integrable. 



11 Discussion 

The lesson that we can learn from this analysis is that the characteristic features of 
the Fermi-Pasta-Ulam chain, such as quasiperiodicity and nonergodicity, are not just 
a property shared by all low-energy solutions of resonant Hamiltonian systems. On 
the contrary: the periodic FPU chain is a rather special system possessing particular 
symmetries and eigenvalues. These cause or may cause nondegenerate integr ability of 
the Birkhoff-Gustavson normal form of the chain, which in turn implies that the KAM 
theorem (cf. is applicable. Still, some questions remain unanswered: 



1. From corollary 9.1 we know that the truncated normal form of the odd FPU 



chain is integrable. In section 10.1 we checked a nondegeneracy condition for the odd 
/3-chain and were able to apply the KAM theorem. Can the truncated normal form of 
the odd chain explicitly be computed also if a ^ 0? Is it really nondegenerate, as we 
are tempted to assume? 

2. What is the reason that the truncated normal form of the even /?-chain is inte- 
grable as we know from corollary |9.3[ ? Is there some hidden symmetry-like property 
of the FPU chain that prevents terms bjbk (j ^ k) from appearing in the truncated 
normal form, thus causing the integrability? 

3. Is it possible to explicitly construct action-angle coordinates for the truncated 
normal form of the even /3-problem, globally or locally, and verify the KAM nondegen- 
eracy condition? 



4. What about the even a-/3 chain? As indicated in corollary 9.3 its truncated normal 
form has a lot of conserved quantities. But is it also really Liouville integrable? If 
yes, then there is a big chance for the KAM theorem to work. And otherwise: can 
we find a counterexample of an even a- (5 chain with many ergodic orbits of low energy? 

Where the second question is of a rather philosophical nature, the other three involve 
tough computations. Answers might be given in a subsequent paper. 
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A Proof of theorem gTt 



This appendix is based on notes of Frits Beukers. Its main intention is to prove theorem 



8.1. Some algebra is used that might be uncommon to the reader, but fortunately the 



conclusions of theorem 8.1 and theorem 8.2 are easily understood 



A.l Sums of roots of unity 

We are interested in solving the resonance equation u(Q, 6) = 0, that is we want to 
find vanishing sums of the eigenvalues ±iuij = ±2isin(^). A study of these sums is 
possible if we first consider sums of roots of unity. 

Fix N £ N. We study the equation Ci + C2 + • • • + (n = in the unknown roots 
of unity Q. The solutions will be determined modulo permutation of the terms and 
multiplication by a common root of unity. We also assume that there are no vanishing 
subsums, that is Yliei d f° r all J C {1, ... , N}, |/| < N. We first state our basic 
tool. Let K be a field generated over Q by roots of unity. Let p k be a prime power 
and let C := e 2ni / pk . Suppose ( K and ( p G K. 

Proposition A.l The minimal polynomial of C over the field K is given by X p — C p 
if k > 2 and X p - 1 + X p ~ 2 + • • • + X + 1 if k = 1. 

For the proof of this proposition we refer to [15|, §60-61. 

To return to our problem let us choose M G N minimal so that (d/Cj) M = 1 f° r a U 
i, j = 1, 2, • • • ,N. Since we can multiply every term of our relation with c^" 1 and put 
Ci = Ci/ Ci we may as well assume that all Ci are M-th roots of unity. Let p k be a primary 
factor of M. Set M' = M/p and write Ci = CiC 1 where Ci € K := Q(e 2ni / M ') and 



n, G {0, 1,2, .. . ,p — 1}. Then, according to proposition A.l, the minimal polynomial 
of C over K is X p - Q p if k > 2 and X p ~ x + X p - 2 + • + X + 1 i£ k = 1. 
We now rewrite our relation in the following form 



E E & s = ( R ) 



s=0 rii=s 

If k > 2 the minimal polynomial of C over K is X p — C p - In particular this means that 
there exist no non-trivial i^-linear relations between 1, C, • • • , C p ■ So the relation (R) 
implies that all coefficients are zero, hence Y2 n =s CiC s = for all s = 0, 1, 2, . . . ,p — 1. 
By the minimal choice of M, at least two of the exponents rii,nj should be different. 
Hence the assumption k > 2 leads automatically to vanishing subsums. 

Let us now assume k = 1. Then the minimal polynomial of C over K is A p_1 + 

X p ~ 2 -| \-X + 1. This means that any isT-linear relation between 1, C, • • • , C p_1 must 

have all of its coefficients equal. Hence, (R) implies that all sums 



<p) 



have the same value a. Since we do not want vanishing subsums we necessarily have 
a 7^ 0. This in its turn implies that each of the summations contains at least one term 
and so p < N. This puts a bound on our search range. 
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A. 2 Explicit computations 



In this section we compute vanishing sums of roots of unity having no vanishing sub- 
sums. It should be noted that the solutions are given modulo permutation of terms 
and multiplication by a common root of unity. 

For each of the specific values of N we shall be considering, we denote by M the 
smallest number such that (d/(j) M = 1 for From the previous section we know 
that M is square free and that p < N for all prime divisors of M. Furthermore, we 
also note that if M divides 6, then it is easy to see that the only possible relations 
without vanishing subsums are 1 — 1 = and 1 + 5 + 5 2 = where 5 = e 2m ^ . So we 
shall assume that there is a prime > 5 dividing M. By N > p > 5 the first interesting 
case to be considered is N = 5. 

N = 5. We have 5\M. Then (P) partitions our sum in precisely five parts, each with 
equal sum. Hence 1 + rj + rj 2 + rf + rf = where rj = e 2m ^ 5 . 

N = 6. Then p < 5, hence 5|M. Then (P) partitions our sum in four parts of length 1 
and one with length 2. Hence we see that — 5 — 5 2 + rj + rj 2 + rf + rj 4 = is the 
solution. 

N = 7. Then p < 7. If 7\M then necessarily, 1 + e + e 2 + e 3 + e 4 + e 5 + e 6 = where 

Suppose 5 is the largest prime dividing M. Then (P) gives a partitioning in 31111 
or 22111. The first gives rise to solutions with zero subsums. The second gives rise 
to the solutions (-<5-5 2 )(l+7/)+?7 2 +?7 3 +?7 4 = and (— 5— 5 2 )(\+rj 2 )+ri+rf' + r q 4: = 
0. 

N = 8. Then p < 7. If 7\M then (P) implies that we have a partitioning 2111111 and 
-5 - 5 2 + e + e 2 + e 3 + e 4 + e 5 + e 6 = 0. 

Suppose 5 is the largest prime dividing M. Then (P) gives a partitioning 41111, 
32111 or 22211. The first two give rise only to vanishing subsums. The last 
solution gives rise to (—5 — 5 2 )(l + rf + rf) + r\ k + rf = where {i,j, k, 1} = 
{1,2,3,4}. 

A. 3 Sums of the iuj 

We are interested in vanishing sums of the eigenvalues ±iujj = ±2isin(^). So we 
look for all solutions of £i + ■ ■ ■ + Cn = such that together with each Q, minus its 
complex conjugate — d" 1 also occurs. Since we shall only be interested in sums of 3 or 
4 eigenvalues iuj, we restrict ourselves to N = 6,8. We include sums with vanishing 
subsums, except vanishing subsums of the form ( — £ = 0, since these give rise to 
vanishing subsums of ioj,fs. So all vanishing subsums of roots of unity must have 
length at least three. 

N = 6. To bring our relation without zero subsums in the desired form, we have to 
multiply it by ±i and we derive 

2isin(vr/6) + 2isin(7r/10) - 2i sin(37r/10) = . 

Now we look at relations with vanishing subsums. There can only be two van- 
ishing subsums of length three. Hence (d + C2XI + 5 + 5 2 ) = with (1,(2 
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arbitrary. It is necessary and sufficient to assume that C1C2 = ~1- This means 
(C — + 5 + S 2 ) where ( is arbitrary. Hence, 

2t sin(vrr) + 2« sin(vr(r + 2/3)) + 2* sin(vr(r + 4/3)) = , 

where r is an arbitrary rational number. 
N = 8. Let us first see what we get from our relations without zero subsums. We find 

2i sin(vr/6) + 2% sin(3vr/14) - 2i sin(vr/14) - 2i sin(5vr/14) = 
2isin(vr/6) + 2i sin(13vr/30) - 2i sin(7vr/30) - 2i sin(3vr/10) = 
2isin(7r/6) + 2% sin(7r/30) - 2% sin(llvr/30) + 2i sin(vr/10) =0 . 

Any relation with vanishing subsums must have subsums both of length 4, or 
subsums of lengths 3 and 5. The first case cannot occur, but the second yields 
+ 5 + 5 2 ) + £2(1 + 77 + • • • + 7/ 4 ) = 0. Both £1, £2 must be purely imaginary 
and have opposite sign. So we can take £1 = — (2 = *, hence 

2% sin(7r/2) - 2i sin(vr/6) + 2i sin(vr/10) - 2% sin(3vr/10) = 



A. 4 Proof of theorem 



We indicate how theorem |8.1| can be proved using the previous paragraphs. 

From the first relation in section IA.3; we infer that icon + icoiL — iuiz-a = if n 

6 10 10 

is a multiple of 30. So multi-indices Q,9 can be found such that |0| + \9\ =3 
and ±z^(0, 6) = iujr + iujn. — icuzn = 0. But for this and 9, we must have that 



yu(0, 6) = ±jr db i °f which one easily verifies that it is unequal to modulo n. 

One finds the same result for the other third order relation of the previous section. 
The verification is not hard, but needs more bookkeeping because of the appearance of 
the arbitrary rational. The conclusion is that for all multi-indices Q,9 with |0| + |0| = 3 
and v{®, 9) = 0, we have that ^(0, 9) 7^ mod n. This proves the first part of theorem 



8.1, which actually states that P3 n ker ad# 2 is too small to have a nontrivial intersec- 



tion with ker(T* — Id). 

The proof of the second part of theorem |8.1| is not harder. For |0| + \9\ = 4, there 
are a number of trivial solutions to the equation ^(0, 9) = 0, namely those of the form 
iujj — iujj + iujk — iuJk = 0. These give rise to the and 9 mentioned in theorem |8.1|. All 



the other solutions to v = must be of the form mentioned in section A. 3 under the item 
l N = 8'. From the first relation we see for instance that ilor + iu3n — iujjL — iivsn = 

6 14 14 14 

if n is a multiple of 42. So multi-indices 0, 9 with |0| + \ 9\ =4 can be found such that 
±z/(0, 9) = iujn + iujsn — iuiiL — iuobn = 0. But for these multi- indices, one must have 

V ' 6 14 14 14 

that ^(0, 9) = ±g ± || ± ± I? mod n. One finds the same conclusion for the 



other relations under the item l N = 8'. This poves the second part of theorem 3.1 
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